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Predictions  of  atmospheric  transmittance  in  desert  aerosol  environments  using  MODTRAN  code  diverge  significantly  from  measur^  data.  Good 
prediction  of  the  desert  particulate  size  distribution  is  required  in  order  to  predict  atmospheric  scattering  and  absorption  parameters.  It  is  also  essential 
to  the  prediction  of  the  aerosol  atmospheric  modulation  transfer  function  which  is  often  the  dominant  component  of  the  overall  atmosphenc  mik 
Recently  an  effort  to  predict  statistics  but  not  size  distribution  according  to  simple  weather  parameters  has  been  made  for  coarse  desert  aerosols.  A 
quantitative  analysis  of  the  desert  particulate  size  distribution  models  was  also  performed.  In  this  research,  the  size  distribution  parameters  measured 
by  optical  counters  are  related  to  weather  parameters.  Known  statistical  and  analytical  models  such  as  MODTRAN  relate  the  size  distnbi^on 
parameters  only  to  relative  humidity  for  continental  atmospheres.  Although  humidity  has  a  significant  role  In  the  prediction  of  aerosol 
other  weather  parameters  are  seen  here  to  strongly  influence  also  the  size  distribution  parameters.  Comparisons  such  as  the  above  can  be  usea  o 
predict  under  which  conditions  the  MODTRAN  aerosol  models  have  good  or  poor  accuracy.  It  is  also  hoped  that  they  will  lead  to 
MODTRAN,  improving  the  accuracy  of  the  humidity  dependence  as  well  as  by  incorporating  other  meteorological  parameters  into  the  MUDTRAN 
prediction  models. 
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Prediction  of  Atmophseric  MTF  and 
Application  to  Image  Restoration,  Based  on  Meteorological  Data: 

First  Quarterly  Report 


This  quarter  work  has  been  concentrated  on  image  restoration  techniques  based 
upon  atmospheric  moduladon  transfer  function  (MTF),  as  well  as  ordering  equipment 
and  getting  organized. 

The  image  restoration  work  is  intended  for  near  infrared  wavelengths  according 
to  the  laser  radar  wavelength  (1.06|im)  at  Eglin  AFB.  At  such  wavelengths  image 
degradation  by  the  atmosphere  derives  from  both  turbulence  and  aerosol  forward 
scatter.  The  former  is  dominant  usually  around  mid-day  while  the  latter  is  usually 
dominant  the  rest  of  the  time. 

A  novel  image  restoration  technique  has  been  developed  here  which  utilizes  a 
Weiner  filter  optimized  for  turbulence  effects,  which  are  often  strongly  time-varying. 
This  technique  corrects  also  for  aerosol-derived  blur,  unlike  adaptive  optics  techniques 
which  correct  only  for  time-varying  tilt  (turbulence).  The  technique  is  described  in  the 
enclosed  appendix.  It  is  very  quick  and,  with  modem  parallel-processing  transputers, 
is  capable  of  producing  restored  images  in  a  fraction  of  a  second.  This  new  technique 
is  based  upon  knowledge  of  average  MTF.  In  the  experiments  here,  average 
atmospheric  MTF  was  measured. 

Efforts  in  the  next  few  quarters  are  to  be  concentrated  upon  prediction  of 
atmospheric  MTF  according  to  weather,  thus  permitting  image  restoration  according  to 
weather.  A  technique  for  predicting  turbulence  MTF  has  already  been  developed. 
Hence,  efforts  are  to  be  concentrated  upon  prediction  of  aerosol  MTF.  To  this  end,  a 
pc  version  of  Modtran  (USAF  Geophysics  Lab)  has  been  purchased  so  as  to  enable  us 
to  work  with  marine  atmospheres  when  winds  at  Eglin  are  from  the  south.  The 
MODTRAN/LOWTRAN  marine  models  are  known  to  be  fairly  reliable.  We  are 
developing  our  own  models  to  predict  aerosol  size  distribution  and  MTF  for  continental 
atmospheres.  These  are  to  be  compared  to  Modtran  predictions  and  to  experimental 
measurements. 

Transmission  predictions  are  to  derive  from  aerosol  size  distribution 
predictions,  as  well  as  aerosol  MTF  predictions. 


APPENDIX 


High  resolution  restoration  of  images  distorted  by  the  atmosphere, 
based  upon  average  atmospheric  MTF 


D.  Sadot,  A.  Rosenfeld,  G.  Shuker,  and  N.  S.  Kopeika 
Department  of  Electrical  and  Computer  Engineering 
Ben-Gurion  University  of  the  Negev 
Beer- She  va,  Israel 


Abstract 

A  new  method  of  real-time  high  resolution  imaging  through  the  atmosphere  is 
presented.  This  technique  is  based  on  the  knowledge  of  average  atmospheric  MTF  at  the  time 
the  image  is  received.  Atmospheric  effects  are  modeled  by  a  noisy  spatial  frequency  filter 
including  an  average  component  described  by  the  average  atmospheric  modulation  transfer 
function,  and  a  noisy  component  modeled  by  the  atmospheric  point  spread  function's  power 
spectral  density.  Analytical  results  are  accompanied  by  experimental  image  restoration 
examples,  indicating  significant  image  quality  improvement  based  upon  knowledge  of  average 
atmospheric  MTF.  This  method  can  be  used  to  help  overcome  the  jitter  characteristics  of 
turbulence,  and  is  capable  of  yielding  real-time  image  restoration  with  resolution  1imitp.fl 
essentially  only  by  the  hardware  itself 


1.  Introduction 


Imaging  through  the  atmosphere  has  progressed  significantly  in  the  last  decade,  in  both 
the  visible  and  infrared  spectral  ranges.  It  is  the  atmosphere  which  usually  limits  image  quality, 
particularly  for  long  atmospheric  paths.  The  main  atmospheric  distortions  are  caused  by  optical 
turbulence,  and  scattering  and  absorption  by  particulates  in  the  atmosphere.  The  main  effect  of 
the  turbulent  medium  over  long  exposures  is  to  produce  wavefront  tilt,  which  causes  image 
shifts  at  the  image  plane.  These  can  be  whole  image  shifts,  or  different  shifts  of  different  parts 
of  the  image,  depending  on  the  isoplanatic  patch,  according  to  the  turbulence  parameters 
(turbulence  strength,  inner  and  outer  scales).  The  image  distortions  caused  by  the  wavefront 
tilt  (typically  on  the  order  of  tens  or  maybe  hundreds  of  microradians)  can  be  partially 
compensated  for  either  by  adaptive  optics  techniques  or  by  using  sufficientUy  short  exposure 
time,  less  than  the  characteristic  fluctuation  time  (usually  a  few  milliseconds).  Another  way  is 
by  using  a  wavefront  sensor  to  construct  the  image  phase  in  addition  to  image  intensity,  and 
therefore  restore  the  image  in  a  deterministic  wayi.  This  is  however  not  practical  for  real  time 
restoration.  Turbulence  distortion  effects  are  characterized  in  botii  short  and  long  exposures 
by  modulation  transfer  functions  (MTFs)2-3.  These  are  statistical  averages  of  turbulence 
random  processes  which  have,  in  addition,  a  nonnegligible  short-time  variance  which  affects 
significantly  image  blur.  In  addition  to  turbulence,  there  are  scattering  and  absorption  effects 
produced  by  molecules  and  aerosols  in  the  atmosphere.  These  cause  both  attenuation  and 
image  blur,  according  to  the  atmospheric  aerosol  MTF^-^.  Unlike  turbulence,  the  aerosol  MTF 
affecting  the  image  is  also  related  strongly  to  instrumentation  limitations^’'^on  angles  of  light 
scatter  actually  recorded  in  the  image.  Thus,  aerosol  effects  can  be  approximated  by  an  MTF 
of  limited  long-time  variation,  according  to  weather  conditions. 


2.  Theory 


Here,  a  new  method  is  developed  for  restoration  of  images  distorted  by  the 
atmosphere,  based  on  prior  knowledge  of  both  average  turbulence  and  aerosol  MTFs.  The 
atmosphere  is  modeled  as  a  random  filter 


h’=h  +  ni  (1) 

where  h'  is  the  instantenous  atmospheric  point  spread  function  (PSF),  h  is  the  average 
atmospheric  PSF,  and  M;  is  an  additive  random  component  with  zero  expectation.  Using  this 
model,  the  image  received  at  the  imaging  system  after  propagating  through  the  atmosphere  is: 

g(x,y)=(/i(x,y)-l-7ii(j:,y))0/(x,y)  +  n2(x,y)  (2) 

where  g  is  the  received  image, /is  the  object,  x  and  y  are  transverse  spatial  coordinates,  and  n2 
is  an  additive  noise  imposed  by  the  instrumentation,  including  optics,  digitization,  electronics 
etc.,  but  not  by  the  atmosphere.  Fourier  transforming  (2)  yields: 

G{u,v)  =[H(u,v)  +  Ni(u,v)]‘  F(u,v)  +  N2iu,v)  (3) 

where  G,  H,  F,  Nj  and  N2  are  Fourier  transforms  of  g,  h,f,nj  and  112  respectively,  and  u  and  v 
are  spatial  frequency  coordinates.  The  received  image  thus  is  a  sum  of  a  deterministic  part  G, 
and  a  random  part  N 


G  =  Gi  +  N  (4) 

where  Gi{u,v)  =  H{u,v)- F{u,v)  (4.a) 

N{u,v)  =  F{u,v)-N\{u,v)-\-N2{u,v).  (4.b) 


The  optimal  image  restoration  filter  in  terms  of  maximal  signal  to  noise  ratio  (SNR)  is  the 
Wiener  filter^-^ 


M{u,v) 


_  \H{u,vf\ _ 

H{u,v)  •  (|H(m,v)2  I+[5„„(m,v)  /  v)]) 


(5) 


where  M  is  the  restoring  filter,  and  .S^care  the  power  spectral  densities  (PSDs)  of/ and  n,  n 
being  the  inverse  Fourier  transform  of  N.  Substituting  (4.b)  into  (5)  and  assuming  statistical 
independence  between/  n and  1x2,  yields 


(6) 


M{u,v)  = 


\H{u,vf\ 


H{u,v)\\H{u,vf  +  5^„2«2  («» v)  /  ( M,  v)]) 


The  restored  image  intended  to  represent  the  original  object  is  then  evaluated  by  multiplying 

the  filter  by  the  received  image  spectrum  in  the  frequency  domain: 

\H{u,vf\ 


X  G{u,v) 
F{u,v)  =  --  -  ' 

\  7  /  rr  /  \ 


(7) 


H(»,v)  +  5„2„2(".v)  /  %(",v)]  ■ 

In  order  to  use  (7),  all  its  terms  must  be  either  calculated  or  estimated.  The  term  h{u,v)  is  the 


average  atmospheric  modulation  transfer  function  (MTF).  Assuming  independence  between 
aerosol  and  turbulence  effects,  this  term  can  be  measured  or  calculated  by  a  multiplication  of 

the  turbulence  MTF  (either  short  or  long  exposure  case)  and  aerosol  MTF.  Turbulence  MTF 
can  be  evaluated  with  the  knowledge  of  standard  meteorological  parameters  using  a 

prediction  modeF®  (verified  independently  by  U.S.  Army  Night  Vision  Laboratory)  or 
IMTURB  or  PROTURB,  and  the  aerosol  MTF^  can  be  evaluated  according  to  knowledge  of 

particle  size  distribution,  which  can  also  be  predicted  via  LOWTRAN,  MODTRAN,  or  other 
models^k  The  term  Sjf{u,v)  can  be  estimated  either  by  using  the  received  image  G{u,v)  or 


by  using  estimation  models  of  the  object’s  PSD,  which  has  been  shown  to  obey  a  fractal  model 
in  the  visible  range,  and  a  Markovian  model  in  the  thermal  range^^.  The  term  S„2„2^ii,v)  is 

assumed  to  be  constant  for  all  spatial  frequencies  since  the  additive  noise  ^2  is  assumed  to  be 

white  noise.  This  assumption  is  commonly  used  and  very  practical,  and  it  has  a  relatively  weak 
effect  on  the  Wiener  filter.  The  term  S„,„,(m,v)  is  very  important,  since  it  includes  the  random 

part  of  the  atmospheric  distortions.  If  this  term  can  be  evaluated,  then  a  Wiener  filter  can  be 
used  for  restoration  of  images  distorted  by  the  atmosphere. 

One  way  of  estimating  is  by  a  direct  measurement.  By  using  the  relation 

=  (8) 
where  IT  is  the  Fourier  transform  of  h',  the  variance  of  H'  can  be  evaluated  by  calculating  both 
terms  of  the  right  hand  side  of  (8).  This  can  be  carried  out  by  measuring  a  series  of 
instantaneous  atmospheric  MTFs,  and  evaluating  the  average  of  both  the  MTF  and  its  square. 
This  is,  however,  not  a  very  practical  way,  particularly  when  real  time  image  restoration  is 


concerned.  It  should  be  pointed  out  that  the  contribution  to  the  random  part  of  the  atmospheric 
MTF  is  due  mainly  to  turbulence  rather  than  aerosols.  Therefore,  in  the  rest  of  the  discussion, 
particularly  in  eqs.  (8)-(19),  for  all  practical  purposes  atmospheric  MTF  and  PSF  refer  to 
turbulence  only  since  aerosol  MTF  displays  much  less  jitter.  However,  in  eqs.  (6)-(7)  H{u,v) 

includes  aerosol  MTF  in  addition  to  turbulence  MTF  since  it  refers  to  average  atmospheric 
MTF. 

Here,  an  analytical  way  is  presented  of  calculating  by  evaluating  both  terms 

of  the  right  hand  side  of  (8).  The  second  term  is  the  square  of  the  turbulence  MTF  which 
can  be  predicted^o  or  measured.  The  first  term  e[h''^]  is  to  be  evaluated  here  analyticaUy. 

The  instantaneous  MTF  of  turbulence  is  given  by^ 

jju,.  _  11  ^^V{%{x,y)  +  X{x-Xfu,y-Xfv)]dxdy 

where  %  =  In  —  is  the  wave  log  amplitude  fluctuation,  A  is  the  wave  amplitude  and  Aq  is  the 

wave  amplitude  in  the  absence  of  turbulence,  /  is  the  imaging  system’s  focal  plane,  and  X  is 
the  radiation  wavelength.  In  a  similar  way,  H'^  is  given  by 

11  ^^Vl%ix,y)  +  xix-Xfu,y-Xfv)  +  x{x-Ax,y-Ay)  + 

JJ  JJ  exp[2x(x,y)  +  2x(x-Ax,y-Ay)]dxdydAxdAy 
+X(x  -Ax-Xfu,y-  Ay  - Xp)]dxdydAxdAy 

11  11  exp[2x(x,y)  +  2x(x-Ax,y-Ay)]dxdydAxdAy' 

Under  the  assumption  of  ergodicy,  the  ensemble  average  MTF  will  be  identical  with  the  long 
exposure  MTF.  and  if  we  assume  that  the  spatial  statistics  of  the  medium  are  Wide  Sense 
Stationary  (WSS),  the  expected  values  are  independent  bfx  and  y  and  can  be  factored  outside 
the  integrals.  The  denominator  in  (10)  is  a  normalizing  factor  which  is  the  zero  frequency 
value,  and  will  be  ignored  for  the  rest  of  this  analysis  for  simplicity.  The  result  is  an  average  of 
the  square  of  the  optical  transfer  function  given  by 


E[H''^{u,v)]  =  \\E[tx^[X{x,y)  +  x{x-yu,y-Xjv)-\- 

+X(jc  -  -  Ay)  +  xU  -  Ax  -  Xfu,y  -  Ay  -  yv)]]d^xd^y. 

To  aid  in  this  calculation,  we  use  the  relations^ 

£{exp(2)}  =  exp(£:{z}+Tvar{z}) 
where  z  is  any  Gaussian  random  variable  (as  x  is),  and 

E{x}  =  -cy,{0) 


(11) 

(12) 

(13) 


where  c  is  the  covariance  function.  The  expectation  and  variance  of  the  variable  at  the 
exponent  are: 


^{X(x,y)  +  x(x-VM.y-Vv)  +  X(x-Ax,y-Ay)  + 
+X(x-Ax-X/“>y-Ay-X/v)}  =4x  =  -  4c^(0) 


^var{x(x,y)  +  x(x-X/M,y-Vv)  +  x(x-Ax,y-Ay) 

+X(x-Ax-X/«,y-Ay-X/i^)}  = 


1 

’(%(x,y)  -X)  +  (X(x-Xfu,y-Xjv)-x)  + 

2' 

lE. 

+(X(x-Ax,y-Ay)-x)  + 

.+(%(x-Ax-X/M,y-Ay-Xyv)-x) 

-  2c^  (0)  +  2cy^  (Xfu,Xfv)  +  2c^  (Ax, Ay)  +  (Ax  +  Xfu,Ay  +  Xfv)  + 
+c^(Ax-Xfu,Ay-Xfv). 


(15) 

Substituting  (14)  and  (15)  into  (12),  and  than  into  (11)  yields: 

e{h''^{u,v)]c<\\  exp[-2cj^(0)+2cj^(X/M,Xyv)  +  2c^(Ax,Ay)  + 

(Ax  +  Xfu,Ay  +  Xfv)  +  +c^  (Ax  -  Xfu,Ay  -  Xp)]dAxdAy  = 

=  exp[4cj^  (0)]/J  exp[-  6c^  (0)  +2c^  (Xfu,Xfv)  +  2c^  (Ax, Ay)  + 

+c^  (Ax  +  ,  Ay  +  X/v)  +  ( Ax  -  X/m  ,  Ay  -  Xfv)]dAxdAy. 

Using  the  relation  between  the  turbulence  MTF  and  the  covariance  function^ 


H(u,v)  oc  exp(-Cjj  (0)  +  Cj^  (Xfu,Xp)) 


(17) 
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and  the  change  of  variables: 


Eq.  (16)  can  be  simplified  to  yield: 


(18) 


e[h'\u,v)}  c  h'^  (m,v) JJ  (u'y)Hiu'+uy+v)H(u'-uy-v)du'dv\  (19) 

oo 

Equation  (19)  determines  the  expected  value  of  the  squared  MTF,  or  in  other  words  the  Point 
Spread  Function's  Power  Spectral  Density,  and  is  the  final  result  of  this  analysis.  This 
evaluation  permits  Wiener  filter  restoration. 

The  advantage  of  this  result  is  that  for  this  restoration  technique  the  only  information 

needed  to  evaluate  (19)  is  the  atmospheric  average  MTF.  One  has  to  measure  or  estimate  the 
average  turbulence  MTF,  then  substitute  it  in  (19)  and  (8)  to  yield  and  apply  it  to 

(7).  As  mentioned  above,  the  value  was  evaluated  only  by  considering  turbulence 

and  not  aerosol  effects.  This  is  not  the  case  at  all  when  dealing  with  the  average  value  of  the 
atmospheric  MTF  in  (7).  Aerosols  are  not  negligible  at  all  in  deriving  the  overall  atmospheric 

MTF,  and  sometimes  are  even  the  dominant  part  of  it^^.  Thus,  aerosol  MTF  must  be  included 
in  H{u,v)  in  (7). 


3.  Experiment 

A  set  of  experiments  was  carried  out  in  the  open  atmosphere  m  order  to  examine  this 
method.  An  imaging  system,  including  a  Pulnix  ccd  camera  model  TM-765,  connected  to  a 
Questar  telescope  of  1600  mm  effective  focal  length,  was  located  inside  the  remote  sensing 
laboratory  of  the  electrical  and  Computer  Engineering.  Dept.  Visible  images  were  recorded 
over  a  horizontal  pathlength  of  6.5  km,  with  an  average  path  elevation  of  25  meters.  No 
wavelength  filters  were  used.  Simultaneously,  atmospheric  MTF  was  measured  via  PSF 
measurements.  This  was  carried  out  using  a  4  mW  HeNe  laser  at  0.6328  jxm  wavelength, 
located  adjacent  to  the  objects  m  the  object  plane.  At  6.5  km  distance,  the  laser  aperture  (<  1 


cm)  served  as  a  point  source.  The  output  signal  from  the  ccd  camera  was  digitized  by  a  Data- 
Translation  frame  grabber  and  sent  to  a  486  personal  computer  for  further  analysis.  Fourier 
transform  of  the  PSF  yielded  the  product  of  the  imaging  system  and  atmospheric  MTFs. 
Division  by  measured  system  MTF  (Fig.  1)  yielded  atmospheric  MTF. 

4.  Results 

Typical  measured  average  and  standard  deviation  atmospheric  MTF  for  a  set  of  100 
images  of  the  laser  aperture  are  presented  in  Fig.  2.  Fig.  3  presents  the  average  squared 
measured  atmospheric  MTF  for  this  set.  In  Fig.  4,  the  calculated  average  squared  atmospheric 
MTF  is  presented,  using  (19).  The  ratio  between  the  measured  and  calculated  squared 
atmospheric  MTF  is  presented  in  Fig.  5.  It  can  be  recognized  that  an  excellent  agreement  is 
achieved  between  theory  and  measurement,  except  at  very  high  spatial  frequencies.  This  can  be 
explained  by  small  inaccuracies  in  the  measurement  of  atmospheric  MTF  at  the  very  high 
frequencies,  where  the  system's  MTF  is  very  poor  and  the  MTF  measurement  is  more  likely  to 
be  affected  by  noise.  Division  by  such  low  MTF  values  implies  large  imaging  errors  in 
atmospheric  MTF  even  for  small  errors  in  imaging  system  MTF. 

Examples  of  restored  images  using  this  new  technique  are  presented  in  Figs.  6  and  7. 
The  restoration  was  carried  out  by  using  (7).  The  term  h{u,v)  used  was  the  measured 
atmospheric  MTF.  The  term  S ff  {u,v)  was  estimated  by  best  fit  to  a  fractal  modeF^  -phe  term 
■S'n2n2(“’^)  was  assumed  to  be  white  noise  and  the  term  was  evaluated  via  (19). 

Restoration  time  was  only  about  2  seconds  per  frame.  This  can  be  shortened  to  a  fraction  of  a 
second,  using  parallel  processing  techniques  already  available.  Therefore  restoration  via  this 
method  can  be  in  real  time.  There  is  a  distinct  improvement  in  fine  details  of  the  images,  even 
though  the  image’s  SNR  is  not  degraded  significantly.  This  is  so  in  spite  of  the  severe  imaging 
conditions  (long  horizontal  distance),  where  turbulence  isoplanatic  patch  was  much  less  than 
image  size,  and  standard  Wiener  filters  failed  completely  in  trying  to  restore  the  image,  as  can 
be  seen  in  Fig.  8  where  a  Wiener  filter  characterized  by  (5)  has  been  used.  In  this  case 


Snniu,v)  refers  to  white  noise  only.  On  the  other  hand,  in  the  restored  images  of  Figs.  6  and  7, 

it  is  generally  quite  possible  to  see  detail  at  6.5  km  almost  as  small  as  limited  bv 

.  .  1 

mstrumentation,  i.e.,  Ax  =  - - »  7cm  where  is  system  bandwidth  (=  50  cycles-mrad'^ 

^^max 

from  Fig.  1).  Examples  are  poles  in  a  fence  in  Fig.  6  and  antenna  bars  in  Fig.  7.  This  means 
that,  essentially,  most  of  atmospheric  blur  is  removed  in  the  restoration  process,  and  the  distant 
object  scene  is  observed  as  if  there  were  no  atmosphere.  Such  restoration  based  on 
atmospheric  MTF  does  not  depend  on  target  shape.  It  is  a  fundamental  image  correction, 
which  can  be  followed  by  other  image  processing  techniques  as  desired.  There  are  cases  where 
turbulence  is  not  so  severe,  and  aerosol  MTF  plays  the  most  important  role  in  determining  the 
atmospheric  MTF,  such  as  in  the  case  of  thermal  imaging^^,  and  the  image  restoration  task  is 
much  easier'^  since  the  atmospheric  MTF  then  is  nearly  deterministic.  Here,  however,  the 
correction  is  for  both  aerosol  and  turbulence  derived  blur. 

Although  the  experiment  here  is  based  upon  measured  atmospheric  MTF  via  PSF,  in 
principle  it  can  be  applied  too  to  predicted  atmospheric  MTF.^>io>ii 

5.  Conclusions 

Restoration  of  visible  images  distorted  by  the  atmosphere  is  presented  based  upon 
atmospheric  MTF.  Such  restoration  is  fundamental  and  yields  image  quality  limited  essentially 
by  instrumentation,  as  if  there  were  no  atmosphere.  The  atmosphere  is  modeled  by  a  noisy 
spatial  frequency  filter  including  an  average  component  described  by  the  MTF,  and  a  noisy 
component  modeled  by  the  PSF's  power  spectral  density.  An  analytical  method  of  deriving  the 
turbulence  PSF's  power  spectral  density  is  developed.  This  is  used  to  optimize  Wiener  filter 
signal-to-noise  ratio.  It  seems  that  when  dealing  with  high  resolution  visible  and  near  infrared 
imaging  through  long  distances,  atmospheric  distortions  need  no  longer  be  treated  as  a 
deterministic  filter  represented  by  average  MTF,  but  can  now  be  modeled  as  an  inherent  part  of 
the  noise  in  the  received  image.  Determining  this  random  part  is  essential  to  achieve  improved 
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restored  images,  and  a  method  of  doing  so  is  presented  here.  This  method  can  be  used  to  help 
overcome  the  jitter  characteristic  of  turbulence,  and  is  capable  of  real  time  image  restoration 
via  parallel  processing  transputers. 
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Figure  Captions 


Fig.  1.  Measured  imaging  system  MTF. 

Fig.  2.  Measured  average  atmospheric  MTF  with  standard  deviation. 

Fig.  3.  Measured  atmospheric  PSF's  power  spectral  density. 

Fig.  4.  Calculated  atmospheric  PSF's  power  spectral  density. 

Fig.  5.  Division  of  measured  atmospheric  PSF's  power  spectral  density  by  calculated 
atmospheric  PSF's  power  spectral  density. 

Fig.  6.  Comparison  between  original  and  restored  atmospherically  distorted  images  for  6.5  km 
distance. 

Fig.  7.  Same  as  in  Fig.  5,  but  for  a  different  scene  at  same  distance. 

Fig.  8.  Image  restoration  using  standard  Wiener  filter. 


System  MTF 


0.3 


0.2 


0.1 


qI - > - 1 - 1 - 1 _ I _ I _  I 

0  5  10  15  20  25  30  3 

Spatial  frequency[cyc!es/mrad] 


i 


measured  PSF  spectral  density 


x10 


spatial  frequency  (1/rad) 


xIO 

spatial  frequency  (1/rad) 


calculated  PSF  spectral  density 


F'J  V 


original 


restored 


I 


restored 


! 

I 


•••:•  ;  .-.Av 


